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Abstract 

By suitably generalizing the Fourier constraint projection in the difference map phas- 
ing algorithm, an object can be reconstructed from its diffraction pattern even when 
the latter has been incoherently averaged over a discrete group of symmetries. This 
resolves an ambiguity in the recent proposal for aligning molecules by means of their 
anisotropic dielectric interaction with an intense light field. The algorithm is demon- 
strated with simulated data in two and three dimensions. 

1. Introduction 

If one is primarily interested in molecules, then crystals are merely a convenient align- 
ment mechanism for improving the signal in a diffraction experiment. The near perfect 
alignment of molecules within a crystal, however, comes at a cost. In a crystal, align- 
ment is achieved by introducing translational periodicity as well, and this restricts the 
sampling of the diffraction pattern to the crystal's reciprocal lattice. A more complete 
sampling of the diffraction pattern, also called oversampling, is desirable when retriev- 
ing the phases in the pattern. This has led to the suggestion of alternate alignment 
mechanisms (Spence & Doak, 2004), where not only would the molecules lack transla- 
tional periodicity, but by arranging the molecular separations to be outside the range 
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of coherence of the illuminating radiation, all inter-molecular interference would be 
eliminated. If such an "incoherent alignment" mechanism were realized, the scattered 
radiation would form the diffraction pattern of an individual molecule and be signal 
enhanced in proportion to the number of molecules in the illuminating beam. 

One of the leading candidates for an incoherent alignment mechanism is the interac- 
tion of an elliptically polarized light field with the anisotropic molecular polarizability 
(Larsen et al. , 2000). As already noted (Spence et al. , 2005), this scheme suffers from 
the flaw that the energetics of the molecular orientation is degenerate with respect 
to a discrete symmetry group of order four generated by 7r-rotations about the prin- 
cipal axes of the polarizability tensor. The scattered radiation from many molecules 
will therefore be an average of four, generally distinct diffraction patterns, all related 
by symmetry. An incoherently averaged pattern will have no simple (linear) relation- 
ship to the scattering density: in particular, it does not correspond to the symmetry 
average of the scattering density. Incoherence proves to be a blessing here as well, how- 
ever. As we show below, the molecule's scattering density can be reconstructed from 
a symmetry-averaged diffraction pattern by a relatively straightforward modification 
of a projection-based phase retrieval algorithm. 

Not surprisingly, reconstruction from a symmetry-averaged diffraction pattern places 
greater demands on the degree of oversampling of the data. This should not pose a 
problem for the proposed alignment mechanism when applied to molecules, however, 
since the oversampling is then not limited by the transverse coherence of the illumi- 
nating beam, but rather, the time needed to acquire the signal. 

2. Difference map algorithm 

We use the difference map algorithm (Elser, 2003a & 2003b), a general purpose itera- 
tive method for finding a point in a Euclidean space that lies in the intersection of two 
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constraint sets. The key components of the algorithm are the two projections, that 
given an arbitrary input, return as outputs the nearest points on the corresponding 
constraint sets. For the reconstructions considered here, the first projection imple- 
ments the support and positivity constraint: 

n ^ , = f if r S or p r < 0, 

s+ Pr Pr j ^ otherwise. 

Here p r is the density at the direct space point r and S is the known or estimated 

support of the scattering density. 

The second projection implements the constraint provided by the diffraction data 

and is most easily expressed in terms of the density in Fourier space, p q . Since we 

have a discrete group G that averages the diffraction pattern, we consider the action 

of this projection on an orbit of densities {p gq } g eG '■ 

n F : pga >-» Pgd ■ (2) 

The constraint is given by 

El/Wl 2 = ^ q , (3) 

where Icq is the diffraction intensity associated with the orbit Gq. The output of IIf 
satisfies (3) while also minimizing the distance 

IIp-pT = E El/v^/v/l 2 - (4) 

Gq 5 eG 

The derivation of a formula for the projection Hp proceeds in two steps. First, 
consider the nonnegative amplitudes 

fgn = \Pgq\ ( 5 ) 

/W = I/WI- (6) 

For any value of / sq ' (to be determined later), the complex Fourier amplitude pgq 
that minimizes (4) subject to (6) is given by the projection to the circle: 

P 9q • (7) 
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The distance (4) can now be expressed in terms of the amplitudes: 

IIp-pT = E E(/, q -0 2 - (8) 

Gq g&G 

In the second step of the derivation we minimize (8) with respect to / 9q ' subject to 
the constraint (3) written as 

E(/W) 2 = ^ q - (9) 

g&G 

This is solved by projecting to the sphere, in analogy to the circle projection in the 
complex plane, only now the dimension of the sphere is one less the cardinality of the 
orbit Gq (a divisor of the order of G): 

^Gq 

{jlgecUg^Y 
Combining this with (7) we obtain the formula 

nH/W = ~, ,i/2 Pa<\ ■ ( n ) 

[EgeG \pg<i\ 2 ) 

Below we use the difference map with (3=1, for which it reduces to a generalization 
of Fienup's input-output map (Fienup, 1982): 

D:p^p' = p + A(p) (12) 

A(p)=U s+ (2U F (p)-p)-U F (p) . (13) 

When the diffraction pattern is not symmetry- averaged, the limiting oversampling 
ratio for real- valued objects is a > 2 (Miao et al. , 1998). With symmetry-averaging, 
the number of Fourier constraints is reduced by the order of the group and the limit 
becomes a > 2\G\. For the alignment scheme considered above, |G| = 4. 



3. Numerical experiments 

The diffraction pattern of a real- valued object in two dimensions has inversion sym- 
metry and results in a two-fold ambiguity in the reconstructed object. When the 
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diffraction pattern is averaged with respect to a mirror, the ambiguity becomes four- 
fold and is generated by two orthogonal mirrors. We demonstrate this for the object 
"R" , one of the least symmetrical letters in the alphabet. Figure 1 shows the object in 
a 27 x 27 pixel support, the diffraction pattern when expanded into a 128 x 128 field 
of view, and the diffraction pattern averaged with respect to a vertical (or horizontal) 
mirror. The algorithm used a slightly enlarged, 30 x 30, support constraint to recover 
the object. Figure 2 shows the time series of the difference map error metric ||A|| given 
by (13). When ||A|| vanishes, the density 

Psol = n 5+ (2n F ( P ) - P ) (14) 

satisfies the support/positivity constraint but also the symmetry-averaged Fourier 
constraint, since 

U s+ (2U F (p)-p) = U F (p) . (15) 

Also shown in Figure 2 is the reconstructed object; in this case the mirror of the 
original. 

Our three dimensional reconstruction experiment reproduced the symmetry am- 
biguity of the proposed molecular alignment mechanism (Spence et al. , 2005). A 
real-valued density was generated by applying a rectifier (positivity projection) to a 
random, low-pass filtered density on a 12 x 12 x 12 support. The diffraction pattern 
was computed for a 32 x 32 x 32 field of view and averaged with respect to the four 
element group generated by 7r-rotations about axes along the sampling grid. When 
combined with the Friedel symmetry of the non-averaged diffraction pattern, the input 
to the reconstruction algorithm thus has an eight element symmetry group generated 
by the three coordinate reflections. 

Figure 3 shows the times series of the difference map error metric for ten successive 
reconstructions from different random starts. The corresponding densities are shown 
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row-wise in Figure 4. Of the eight possible symmetry related reconstructions, the ten 
trials produced six. 

Acknowledgement 

This work was supported by Department of Energy grant DE-FG02-05ER46198. 

References 

Elscr, V. (2003a). J. Opt. Soc. Am. A20, 40-55. 
Elscr, V. (2003b). Acta Cryst. A59, 201-209. 
Fienup, J. R. (1982). Appl. Opt. 21, 2758-2769. 

Larsen, J. J., Hald, K., Bjerre, N., Stapelfeldt, H. & Seideman, T. (2000), Phys. Rev. Lett. 85, 
2470-2473. 

Miao, J., Sayre, D. & Chapman, H. N. (1998). J. Opt. Soc. Am. A15, 1662-1669. 

Spcnce, J. C. H. & Doak, R. B. (2004). Phys. Rev. Lett. 92, 198102-198104. 

Spcncc, J. C. H., Schmidt, K., Wu, J. S., Hcmbrcc, G., Weierstall, U., Doak, B. & Fromme, 
P. (2005). Acta Cryst. A61, 237-245. 



IUCr macros version 2.0/35: 2001/06/20 



7 




Fig. 1. Left: object of study in a 2D reconstruction experiment. Center: diffraction 
pattern, right: diffraction pattern averaged with its mirror. 




Fig. 2. Left: Time series of the difference map error metric ||A|| vs. iteration number; 
right: reconstructed object. 
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Fig. 3. Error metric time series for the ten 3D reconstructions shown in Figure 4. All 
solutions required fewer than 300 iterations. 
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Fig. 4. Each horizontal row of twelve, 12 x 12 pixel arrays, represents one 3D recon- 
struction from symmetry averaged data. Of the eight symmetry related reconstruc- 
tions possible, the ten experiments shown realize six of these (rows resorted to ease 
comparison). 



Synopsis 

An object can be reconstructed from its diffraction pattern even if the latter is incoherently 
averaged with respect to a discrete group of symmetries. 
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